arXiv:1503.05728v2 [hep-ph] 16 Sep 2015 


Higher radial and orbital excitations in the charmed meson family 


Qin-Tao Song^’“’"^,* Dian-Yong Chen*’^^,* Xiang and Takayuki Matsuki^’®” 

'Nuclear Theory Group, Institute of Modern Physics of CAS, Lanzhou 730000, China 
^Research Center for Hadron and CSR Physics, Lanzhou University & Institute of Modern Physics of CAS, Lanzhou 730000, China 
^School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China 
''University of Chinese Academy of Sciences, Beijing 100049, China 
^Tokyo Kasei University, I-18-1 Kaga, Itabashi, Tokyo 173-8602, Japan 
^Theoretical Research Division, Nishina Center, RIKEN, Saitama 351-0198, Japan 

Having abundant experimental information of charmed mesons together with the present research status, we 
systematically study higher radial and orbital excitations in the charmed meson family by analyzing the mass 
spectrum and by calculating their OZI-allowed two-body decay behaviors. This phenomenological analysis 
reveals underlying properties of the newly observed charmed states £1(2550), £>*(2600), £>*(2760), £>(2750), 

£>j(2580), £)}(2650), £>j(2760), £)j(2740), £)j(3000) and £>}(3000) to provide valuable information of the 
charmed mesons still missing in experiments. 


PACS numbers: 14.40.Lb, 12.38.Lg, 13.25.Ft 

I. INTRODUCTION 

In the past years, experimentalists have made great ad¬ 
vances for observing charmed mesons and accordingly a list 
of charmed mesons collected in Particle Data Group (PDG) 
has become increasingly abundant [ 1 ]. Candidates of higher 
radial and orbital excitations in the charmed meson family are 
the newly observed mesons D{2550), D*(2600) D*(2760) and 
D(2750) from BaBar [2], and Dy(2580), D}(2650), D}(2760), 
Dy(2740), Dy(3000) and D}(3000) from LHCb [3], whose ex¬ 
perimental information together with other observed charmed 
mesons is given in Table I. 

This new observation and the present research status mo¬ 
tivate us with great interest to carry out systematic and phe¬ 
nomenological study of higher radial and orbital excitations 
in the charmed meson family, which will reveal the underly¬ 
ing properties of the observed charmed mesons and provide 
much more information for their further experimental search. 

Although there have been a couple of works studying the 
heavy-light systems including charmed mesons with their de¬ 
cay modes [4-11], in this work we focus mainly on applica¬ 
tion of our modified Godfrey-Isgur model developed in Ref. 
[12] to calculating mass spectrum and decay behaviors of 
charmed mesons. Difference of our model from the Godfrey- 
Isgur model (GI) proposed in Ref. [13] is that the confining 
potential br in the GI model is replaced by the screened one 
to reflect the unquenched effect (see Ref. 
[12] for more details). This model includes such a screening 
effect that a linear confinement term br be screened or soft¬ 
ened at large distances by virtual quark pairs and dynamical 
fermions [14, 15], which is important especially for higher ra- 
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dial and orbital excitations. Even though the GI model does 
not consider the screening effect, it has been successful to de¬ 
scribe the low-lying states [13]. With the modified GI model, 
we will show a quite successful mass spectrum of the high- 
lying charmed meson family. Having the abundant experi¬ 
mental data, we are able to compare our theoretical values 
with the corresponding experimental results, which can not 
only test the reliability of our modified GI model, but also give 
some useful structure information of the observed charmed 
mesons. In addition, we will predict masses of some missing 
charmed mesons, which gives an important hint to experimen¬ 
tally explore these missing states. 

One more valuable lesson we learned from the former work 
Ref. [12] is that an experimental value of a mixing angle 
which is close to the one in the heavy quark limit determines 
relative decay widths, broad or narrow, with the same 
states in two spin multiplets. These results are found to be 
model-independent. The same scenario can be applied to the 
charmed meson family too and details will be described later. 

After analyzing the mass spectrum of charmed mesons, in 
this work we further study two-body OZI-allowed decays of 
charmed mesons, where the quark pair creation (QPC) model 
[16-22] is applied to calculating strong decays. With the mod¬ 
ified GI model, we obtain the numerical wave functions of 
charmed mesons, which can be applied to calculating strong 
decays of charmed mesons, where we give partial and total 
widths of charmed mesons under discussion and some typical 
ratios relevant to these decays. 

Combining the analysis of mass spectrum with calculation 
of the decay widths, we can learn the properties of higher ra¬ 
dial and orbital excitations in the charmed meson family. Fur¬ 
thermore, we can also shed light on the underlying structures 
of the observed charmed states under discussion, which is the 
main task of this work. 

This paper is organized as follows. In Sect. II we give a 
concise review on the observed charmed mesons. Then, the 
mass spectrum will be analyzed in Sect. Ill, where the modi¬ 
fied GI model is briefly introduced. In Sect. IV we will give 
a brief review of the QPC model and calculate the two-body 
OZI-allowed decay of the charmed mesons under discussion 
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TABLE I: Experimental information of the observed charmed mesons. 


State 

Mass (MeV) [1] 

Width (MeV) [1] 

T" observation 

Observed decay modes 

D 

1864.84 + 0.07 




D- 

2010.26 + 0.07 

0.083 



Di(2420) 

2421.4 + 0.6 

70 + 21 

ARGUS [23] 

D*7T 

D;(2460) 

2464.3 ± 1.6 

20+10 + 5 

TPS [24] 

D*7T^ 

Oi(2430) 

2427 ±26 + 25 

384+*“’ ± 75 

Belle [32] 

D'n 


2477 ± 28 

266 ± 97 

BaBar [38] 

D*7T 

0^(2400) 

2308 ±17 + 32 

276 + 21 +63 

Belle [32] 

Dtt 


2407 ±21 + 35 

240 ±55 + 59 

EOCUS [39] 

Dtt 


2297 ±8 + 20 

273 ±12 + 48 

BaBar [40] 

Dtt 

D(2550) 

2539.4 ± 4.5 ± 6.8 

130+ 12+ 13 

BaBar [2] 

D’tt 

D*(2600) 

2608.7 ± 2.4 ± 2.5 

93 + 6+13 

BaBar [2] 

D**’;r 

D(2750) 

2752.4+ 1.7 + 2.7 

71 + 6+11 

BaBar [2] 

D*n 

£>•(2760) 

2763.3 ± 2.3 ± 2.3 

60.9 + 5.1+3.6 

BaBar [2] 

Dtt 

D(2580) 

2579.5 ± 3.4 ± 5.5 

177.5 + 17.8 + 46.0 

LHCb [3] 

D*7T 

£>•(2650) 

2649.2 ± 3.5 ± 3.5 

140.2+ 17.1 ± 18.6 

LHCb [3] 

D*7T 

£)j(2740) 2737.0 + 3.5 ±11.2 

73.2+ 13.4 + 25.0 

LHCb [3] 

D*7T 

£>•(2760) 

2761.1 +5.1 + 6.5 

74.4 ± 3.4 ± 37.0 

LHCb [3] 

D^'’''tt 

£)y(3000) 

2971.8 + 8.7 

188.1 +44.8 

LHCb [3] 

D*7r 

£);(3000) 

3008.1 +4.0 

110.5 ± 11.5 

LHCb [3] 

Dtt 


by the QPC model. This paper will end with a short summary 
in Sect. V. 


II. STATUS OF THE OBSERVED CHARMED MESONS 

Before investigating the observed charmed mesons, we 
need to briefly review the present research status on candi¬ 
dates of higher radial and orbital excitations of the charmed 
meson family, 


A. £>*(2400), O, (2430), Di (2420) and D*(2460) 

As the first observed P-wave charmed meson, Di(2420) 
was reported in the D*7i invariant mass distribution by the 
ARGUS Collaboration, where its measured mass and width 
are M = 2420 + 6 MeV and T = 70 + 21 MeV, respectively 
[23]. In 1989, the TPS Collaboration confirmed Z)i(2420) in 
the D*^n^ decay channel [24]. According to this decay mode, 
its spin-parity is either 7^ = l^or2^. The angular momentum 
analysis by the ARGUS Collaboration further shows that the 
observed Di(2420) has 7^ = 1+ [25]. The Di(2420) has been 
confirmed by other experiments [2, 26-37] too. 

The Belle Collaboration observed a broad state Di(2430) 
with mass M - 2427 + 26 ± 20 + 15MeV and width T = 
384];'°^ + 24 + 70 MeV by analyzing the B —> D’tttt process, 
where its spin-parity is determined as 7^ = 1 ^ by the helicity 
distributions [32]. In 2006, the BaBar Collaboration studied 
the D*n invariant mass spectrum, where the broad Di(2430) 


was confirmed [38]. 

Besides Di(2430), Belle announced the observation of an¬ 
other broad state Dq( 2400), which has mass M = 2308 + 
17 + 15 + 28 MeV with 7^ = 0+ [32]. The D;(2400) was 
also confirmed by the FOCUS and BaBar Collaborations with 
M = 2403 + 14 + 35MeV [39] and M = 2297 + 8 + 5 + 19MeV 
[40], respectively. Since different experiments gave quite dif¬ 
ferent mass values for Dq( 2400), we will discuss mass depen¬ 
dence of the Dp(2400) decay by varying mass in the range 
(2290 ~ 2350) MeV. The Dg(2400) as a P-wave state with 
jP - 0^ is supported by the theoretical works [4, 5, 8, 41]. 

The TPS Collaboration observed the charmed meson with 
M - 2459 + 3MeV and T = 20+10 + 5 in the invariant 
mass spectrum of [24], which shows that this meson 

has either 7^ = 0^ or 2^. This observation was confirmed 
by the ARGUS Collaboration in the channel, and their 
angular momentum analysis suggests the 7^ = 2^ assignment 
to this state [26]. Thus, this resonance is named ^2(2460). 
Later, the CLEG Collaboration again confirmed the existence 
of ^2(2460), where it has mass M = 2461 + 3 + IMeV and 
width r = 20 ];jq^® 2 MeV. In addition, the ratio 

S(D!(2460) ^ D+7T-) 

- - - =2.3+ 0.8 (1) 

S(D;(2460) ^ D*+7r-) ^ ^ 

was measured as well in Ref. [27], which is consistent with 
the theoretical calculations in Refs. [4, 5, 10, 41-44]. In re¬ 
cent years, many other experiments have reported ^2(2460) 
[2, 25, 28-31, 33-37, 39, 40, 45]. 

The above experimental information indicates that there are 
two charmed mesons with 7^ = 1^. In the heavy quark limit 
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mQ °°,Je ^ L is a good quantum number, where is 
the spin of a light quark and L is its angular momentum. Thus, 
heavy-light mesons can be classified by Two 15-wave 
charmed mesons form one doublet ( 0 ^, 1 “) with = 5 ■ 
Four IP-wave charmed mesons can be grouped into two dou¬ 
blets, ( 0 ^, 1 ^) and ( 1 ^, 2 ^), with and - |^, respec¬ 

tively. Since states with - 1^ in (0^, 1^) and (1^,2^) dou¬ 
blets decay into D*n via S -wave and D-wave [4, 5, 8 , 41, 46- 
48] in the heavy quark limit, respectively, the charmed mesons 
with 7^ = 1 + in ( 0 ^, l"^) and ( 1 ^, 2 +) doublets have broad 
and narrow widths, respectively, which makes us easily distin¬ 
guish two 1 ^ charmed mesons experimentally observed, i.e., 
Di(2430) and Di(2420) belong to the doublets (0^, 1^) and 
(1^, 2^) [4, 5, 8 , 41, 42, 44], respectively. 


B. Z)(2550), Z)j(2580), P*(2600) and £)}(2650) 


The D(2550) was observed by the BaBar Collaboration, 
which has the mass M - 2539.4 + 4.5 + 6.8 MeV and width 
r = 130+12+13 MeV [2]. The D(2550) is suggested to be a 
candidate for D(2'5 0 ) by the helicity distribution analysis [ 2 ]. 
In 2013, an unnatural state 0/(2580) was found in the D*n 
invariant mass spectrum by the LHCb Collaboration through 
the process pp —» DnX [3]. Since the resonance parame¬ 
ters of 0/(2580) are similar to those of 0(2550), 0(2550) and 
0/(2580) are regarded as the same state. 

The mass of 0(2550) or 0/(2580) is consistent with the 
theoretical prediction of 0(2'5 0 ) in Ref. [13]. In addition, 
the decay width of 0(2'5 0 ) was calculated by the QPC model 
[46], which is close to the lower limit of experimental width 
of 0(2550)/0/(2580). In addition, the theoretical studies 
presented in Refs. [49, 50] also show that 0(2550) can be 
O(2'5o), the first radial excitation of O mesons. However, 
the authors in Ref. [5, 51-53] indicated that the theoretical 
total width of 0 ( 2 '5 0 ) is far below the experimental value of 
0(2550). 

The BaBar Collaboration reported a resonance 0*(2600) 
with mass M - 2608.7 + 2.4 + 2.5 MeV and width T = 
93 + 6+13 MeV, which is regarded as a radial excitation of O* 
by the helicity distribution analysis [ 2 ] and they also measured 
the ratio [ 2 ] 


S(0*°(2600) ^ O+TT-) 
S(O*0(2600) ^ 0*+7r-) 


0.32 + 0.02 + 0.09. 


( 2 ) 


A natural state O}(2650) was found in the 0 * 7 r invariant 
mass spectrum by the LHCb collaboration through the process 
pp D*nX [3], where 0^(2650) is tentatively identified as a 
— \^ state, a radial excitation of O*. Therefore, Dy(2650) 
[3] and D*(2600) [2] are the same state. 

In Ref. [54], they predicted that mass of D(2^S 1 ) is 2620 
MeV via the constituent quark model, which is in good agree¬ 
ment with the experimental value of £>*(2600). Moreover, the 
ratio r(D(235i)° ^ n-) IY(D{2^ S xf D’+tt') = 0.47 

was predicted via the relativistic chiral quark model [4], which 
is close to the upper bound of Eq. (2). In Ref.[5], the authors 
calculated a mass spectrum and wave functions of charmed 


mesons via a relativistic quark model, and then adopted the 
obtained masses and wave functions as an input to estimate 
hadronic decay widths. Here, the predicted mass of D(2^S 1 ) 
is 2692 MeV which is heavier than D*(2600), while the pre¬ 
dicted total width of D{2^S 1 ) is consistent with the experi¬ 
mental data of D*(2600). Furthermore, the assignment 2^5 1 
to D*(2600) is also supported by the studies in Refs. [49- 
53, 55, 56]. 


C. Z?*(2760), O)(2760), £)(2750) and Z)/(2740) 


The D*(2760) observed by the BaBar Collaboration in the 
Dn invariant mass spectrum [2] can be assigned to a D-wave 
charmed meson since its mass is consistent with the theo¬ 
retical prediction in Ref. [13]. Later, LHCb announced 
the observation of a natural state Dj{276Q) with mass M = 
2761.1 + 5.1 + 6.5 MeV and width T = 74.4 + 3.4 + 37.0 
MeV. Both D*(2760) and Dy(2760) can be regarded as the 
same state since these two states have the similar widths and 
masses [57]. 

Comparison between the prediction in the Ref. [5] and 
the experimental data of D*(2760) shows that D*(2760) can 
be either D(l^Di) or However, the assignment of 

D*(2760) to D(l^Di) or cannot be supported by the 

result shown in Ref. [46] since the calculated total widths 
of these two assignments are far larger than the experimental 
value. Later, in Ref. [51], it was suggested that D*(2760) is a 
mixture of the 2^5 1 and PDi states, which is also supported 
by the study presented in Ref. [49]. However, calculation 
by the constituent quark model shows that the as¬ 

signment of D*(2760) cannot be excluded [52], which is also 
supported by the works in Refs. [51, 53, 55, 58]. 

Besides D*(2760), another state D(2750) was also observed 
by the BaBar Collaboration in the D* 7 r mass spectrum, where 
its mass and width are M = 2752.4 + 1.7 + 2.7 MeV and 
r = 71+6+11 MeV, respectively [2]. AlthoughD(2750)can 
be a good candidate of a D-wave charmed meson according 
to the mass spectrum analysis in Ref. [13], the helicity distri¬ 
bution analysis of D(2750) does not support the Z)(l^Di) and 
0(1^0^) assignments [2]. BaBar also gave the ratio [2] 


S{D*{21bQf D+tt^) 
S(D(2750)0 ^ D*+;r-) 


0.42 + 0.05 + 0.11. 


(3) 


As an unnatural state, D /(2740) was found by the LHCb 
Collaboration, which has mass M = 2737.0+3.5+ 11.2 MeV 
and - 2“ [3]. Due to the similarity between D(2750) 
and D/(2740), it is possible that D(2750) and D/(2740) are 
the same state. Before the observation of D(2750)/Z)/(2740), 
the masses of two 2“ charmed mesons were predicted in Ref. 
[5], where the masses of the 2“ charmed mesons belonging 
to the (r, 2-) and (2', 3^) doublets are 2883 MeV and 2775 
MeV, respectively, which shows that D(2750)/Z)/(2740) as a 
2“ state with = 5/2“, i.e., (2“, 3“), is more favorable. How¬ 
ever, we also notice that the corresponding theoretical width 
of this 2 “ state is not consistent with the experimental value. 
After observating D(2750)/D/(2740), the authors of Ref. [58] 
calculated the ratio r(D*(2760)° ^ D+;r-)/r(D(2750)° ^ 
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D*^7t~) by adopting an effective Lagrangian approach, which 
is consistent with the experimental data, where £>*(2760) and 
D(2750) are identified as the l^Dj, and 2“ states in the (2 ,3 ) 
doublet, respectively. The studies in Refs. [49, 52, 53, 55] also 
suggested that D(2750) is a iT state in the (2“, 3“) doublet. 


D. £)j(3000) and £)}(3000) 

The LHCb Collaboration observed the unnatural state 
Dy(3000) in the D*ti invariant mass spectrum [3] , where its 
resonance parameters are 

M ^291l.%±%.lUsW, r= 188.1 +44.8MeV. 

Then, different theoretical groups carried out the study of 
Dj(3000). In Refs. [59], Dj(3000) is regard as the first radial 
excitation of Di(2430), which was also confirmed by Refs. 
[56, 60]. However, other possible assignments to Dy(3000) 
were proposed, i.e., the D(3'5'o) [50] and D(3^) [56] assign¬ 
ments. 

A natural state Dy(3000) was also reported by LHCb in the 
Dtt invariant mass spectrum [3], which has 

M = 3008.1 +4.0 MeV, r= 110.5 + 11.5 MeV. 

The D*(3000) is variously explained as D(2^Pq) [59], 
D(I^F 2 ) [56] and D{1^F4) [50, 56, 60]. 


III. MASS SPECTRUM 


For the heavy-light meson system, we need to adopt a rela¬ 
tivistic quark model to study their mass spectrum since a rel¬ 
ativistic effect for a heavy-light meson system is significant. 
The Godfrey-Isgur (GI) model can well describe the meson 
spectrum [13], which is a typical quenched quark model. Af¬ 
ter the discovery of D,o(2317) [61-64], D,i(2460) [62-65] 
and A(3872) [66], theorists realized that it is necessary to take 
into account coupled channel effects, especially for higher ra¬ 
dial and orbital excitations of hadrons [67-70], where the cou¬ 
pled channel effects may change the meson spectrum. This 
motivates us to modify the GI model by considering the cou¬ 
pled channel effects. 

In general, spontaneous creation of light quark-antiquark 
pairs inside a meson can soften a linear confinement poten¬ 
tial br by screening a color charge at distances larger than 
about one fermi [15], which is known as the screening effect. 
The screening effect has been seen by the unquenched Lattice 
QCD and holographic models [71-73]. The mass suppression 
can be caused by both the screening and coupled channel ef¬ 
fects. Although the screening effect can be almost equivalent 
description of the coupled channel effect, we need to empha¬ 
size that the screening effect cannot depict the near-threshold 
effect as the coupled channel effect does [74], which is ob¬ 
tained by studying charmonium spectrum. The authors of 
this paper made a comparison of the results by adopting the 
screening and coupled channel effects. Applying the idea in 


Ref. [74] to our case, we can similarly study a charmed me¬ 
son spectrum. However, it is a complicated task and can be 
assessed at a future theoretical work'. 

In Refs. [12, 75, 76], the screening effect was taken into ac¬ 
count when studying the mass spectra of light mesons, char- 
monia and charmed-strange mesons. Mezzoir et al. provide 
description of a highly excited light-quark meson spectrum by 
flattening the confining potential br at distances larger than r, 
[75]. The screened potential model [77, 78] was adopted to 
compute the charmonium spectrum [76]. In our recent work 
[12], the screening effect was introduced to modify the GI 
model, where the mass spectrum of charmed-strange meson 
family with this treatment is greatly improved compared with 
the results of the GI model. As a sister work of Ref. [12], 
the present work focuses on the charmed mesons applying the 
modified GI model [12] to obtain their mass spectrum. 

In the GI model, the confining potential br is smeared out 
to include relativistic effects, i.e., 

V(r) - cPr'puir — r')br'. (4) 

Here, the GI model is a quenched quark model in the sense 
that the effect of quark-antiquark pairs is not introduced [13]. 
In order to take account of the screening effect in the GI 
model, the confining potential br is replaced with [77, 78] 


br 




b{l-e-n 


(5) 


where y*“(r) behaves like br at short distances and constant 
bjp at large distances. Furthermore, the smearing function is 
introduced to take into account nonlocality property of poten¬ 
tials [13], i.e.. 


= r cPr'pxiir - r')— — -(6) 

J P 

The detailed explanation of how to introduce the screening 
effect into the GI model can be found in Ref. [12]. 

The mass spectrum of charmed mesons is calculated in the 
GI model and the modified GI model, which is shown in Table 
II. Here, all the parameters of the modified GI model are the 
same as those of the GI model [13] except for the additional 
parameter p in y®‘^’^(r). In this Table, we present the theoret¬ 
ical predictions of the mass spectrum with several values of 
p, 0.01 GeV, 0.02 GeV, 0.03 GeV, and 0.04 GeV. In order to 
choose a suitable value for p, we calculate thex^ of predicted 
masses for each p. In Table III we list the theoretical results 
with the least value of p - 0.03 GeV together with experi¬ 
mental data. Here the concrete expression for;^f^ is, 

2 /-^ThCO -'^Exp(0'\^ 


^ We would like to thank the referee for his/her useful suggestion on this 
point 
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FIG. 1: (color online). Mass spectrum of charmed mesons (in units of MeV). Here, the blue lines stand for the results of the GI model [13], 
while the red lines are those of the modified GI model with yu = 0.03 GeV. The purple squares denote the experimental data taken from PDG 
[1], while the green circles and blue lozenges are experimental masses of Ref. [2] and Ref. [3], respectively. The corresponding 
quantum numbers are listed on the abscissa. In addition, when there exists a mixture of and states, we use a notation L^. 


where ^ThiO and ^ExpiO are theoretical and experimental 
values, respectively. Error(/) denotes the experimental error of 
the mass of a charmed meson. A minimum of is obtained 
by our model with fi — 0.03 GeV. Eventually, the modified 
GI model with this value of ^ improves the whole description 
of the charmed meson mass spectrum compared with the GI 
model (see Table III for more details). 

In Table III, we further make a comparison between exper¬ 
imental data and theoretical values obtained via the modified 
GI model. Besides the 2'5'o and 2^S i states, a difference be¬ 
tween theoretical and experimental values of other charmed 
mesons is less than 20 MeV. We can conclude from the mass 
spectrum analysis, i.e., 

1. Two 15 states and four IP states in the charmed me¬ 
son family can be well reproduced by the modified GI 
model. 

2. Both D(2550) reported by the BaBar Collaboration [2] 
and D y(2580) from the LHCb Collaboration [3] are usu¬ 
ally considered as a candidate of D(2'5o). Here, the 
mass of D(2550) is quite close to the theoretical mass 
of D(2'5o), while the mass of Dy(2580) is larger than 
the theoretical mass of D(2'5o) by about 40 MeV. 

3. There exist two experimental results D*(2600) and 
D}(2650) from BaBar [2] and LHCb [3], both of which 


can be a candidate of D(2^S i). Our result shows that 
D*(2600) is closer to the theoretical mass of D(2^S i). 

4. D(2750)/Dy(2740) is a good candidate of D{\D 2 ), 
while D*(2760)/D}(2760) corresponds to Dil^Di) or 
D(PD3). 

5. Dy(3000) and Dy(3000) can be candidates of D(3*5o) 
and D{3^S i), respectively. In addition, the mass spec¬ 
trum study cannot exclude a possibility of D y(3000) and 
D*j{3QQQ) as the IP states in the charmed meson fam¬ 
ily. We also notice that the theoretical masses of the 2P 
states are about 100 MeV smaller than the correspond¬ 
ing experimental data of Dy(3000) and D*j{3000). 

We also notice a study of the heavy-light meson spec¬ 
troscopy within the framework of the QCD-motivated rel¬ 
ativistic quark model based on the quasipotential approach 
[79]. The conclusion of mass spectrum obtained above in the 
present work is consistent with those of Ref. [79]. 

In order to clearly identify the properties of the observed 
charmed mesons, we need to perform a systematic study of 
their decay behaviors, which is the main task in the next sec¬ 
tion. 
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TABLE II: The calculated masses of charmed mesons by the modi¬ 
fied GI model and comparison with those obtained by the GI model. 
Here, we take several /r values, yu = 0.01,0.02,0.03 , and 0.04 GeV 
to show yu dependence of the modified GI model. Values in brack¬ 
ets for the GI model and yu = 0.03 are those of R = 1/yS, which can 
be determined by solving J = J 9(p)^p^d^p, where 

'^nLAf^P^ is an SHO wave function and ibCp) is the wave function of a 
charmed meson which we obtain by solving an eigenvalue equaiton. 
We need to emphasize that we do not consider mixing among states 
with the same quantum number when presenting the results, p is in 
units of GeV, while R is in units of GeV“*. 


GI model 

p = 0.01 

Modified GI model 

p = 0.02 p = 0.03 

p = 0.04 

l‘5o 1874(1.52) 

1869 

1865 

1861(1.54) 

1855 

2'5o 2583(2.08) 

2566 

2550 

2534(2.22) 

2518 

3‘5o 3068(2.33) 

3037 

3005 

2976(2.50) 

2945 

C5i 2038(1.85) 

2032 

2027 

2020(1.89) 

2015 

2^51 2645(2.17) 

2628 

2610 

2593(2.33) 

2576 

3^51 3111(2.38) 

3079 

3047 

3015(2.56) 

2983 

PPi 2457(2.00) 

2447 

2436 

2426(2.08) 

2415 

2'Pi 2933(2.27) 

2909 

2885 

2861(2.44) 

2837 

CPo 2398(1.85) 

2387 

2376 

2365(1.92) 

2354 

2^Po 2932(2.22) 

2907 

2881 

2856(2.38) 

2831 

CP| 2465(2.00) 

2453 

2441 

2431(2.08) 

2419 

2^Pi 2952(2.27) 

2927 

2902 

2877(2.44) 

2852 

CP 2 2501(2.22) 

2490 

2479 

2468(2.33) 

2456 

2^2 2957(2.38) 

2933 

2908 

2884(2.56) 

2859 

1‘D2 2827(2.27) 

2807 

2791 

2773(2.38) 

2755 

2*02 3225(2.44) 

3193 

3160 

3128(2.63) 

3095 

Ed, 2816(2.13) 

2798 

2780 

2762(2.27) 

2744 

2^Di 3231(2.33) 

3198 

3164 

3131(2.56) 

3097 

ED 2 2834(2.27) 

2816 

2797 

2779(2.38) 

2761 

2^D2 3235(2.44) 

3202 

3169 

3136(2.63) 

3102 

EDj 2833(2.38) 

2815 

2797 

2779(2.56) 

2761 

2^D3 3226(2.50) 

3194 

3162 

3129(2.70) 

3097 

I'Fi 3123(2.44) 

3097 

3072 

3046(2.63) 

3019 

EC 2 3132(2.33) 

3106 

3080 

3053(2.50) 

3027 

ETj 3129(2.44) 

3104 

3078 

3051(2.63) 

3025 

EC 4 3113(2.50) 

3088 

3063 

3037(2.70) 

3011 


IV. STRONG DECAY BEHAVIORS 

Study of the two-body Okubo-Zweig-Iizuka (OZI) allowed 
strong decay behaviors of the observed and the predicted 
charmed mesons can provide much more information on the 
features of the charmed mesons under discussion, which in¬ 
cludes total and partial decay widths. What is more important 
is that the inner structure of the observed charmed mesons 
given by the mass spectrum analysis in Sect. Ill can be further 
tested here. 

As an effective approach to study the OZI allowed strong 
decays of hadrons, the quark pair creation (QPC) model is ap¬ 


plied to compute the OZI-allowed strong decays of charmed 
mesons, which was first proposed by Micu [16] and was fur¬ 
ther developed by the Orsay group [ 1 7-22] . Here, a meson de¬ 
cay occurs through a flavor and color singlet quark-antiquark 
pair created from the vacuum. To depict a qq pair creation 
from the vacuum, the transition operator T is introduced as. 


T - -Sy ^(Im; 1 - m|00) j ^fp3^fp4(5^(p3 -i- P 4 ) 

where 7 is a dimensionless constant which reflects the strength 
of creating a quark-antiquark pair from the vacuum. In Ref. 
[80], 7 = 8.7 is obtained for the uTildd pair creation by fit¬ 
ting with the experimental data, where the S U{3) flavor sin¬ 
glet wave function was adopted, while for reflecting the S U(3) 
breaking effect we take 7 = 8.7/ V3 for the ss pair creation as 
suggested in Ref. [17]. We need to explain why there exists an 
extra 1 / V3 factor for the ss pair creation. In Ref. [17], the au¬ 
thors realized the S U(3) and S t/(4) breaking by defining the 
flavor function cpo of a qq pair created from the vacuum, i.e., 
00 reads 0o = mm + dd + crs(_ss) + cr^Ccc) where ctj = 
o'c = niulmc, and cr, expresses the S U{3) breaking. In oder 
to check whether the estimate of ctj is in agreement with 
hadron spectroscopy, the authors of Ref. [17] calculated ra¬ 
tios ra,(i 320 ) = r(a2(1320) ^ 7r;7)/r(a9(1320) ^ KK) and 
runio) = r(/2(1270) ^ 7r7r)/r(/2(1270) ^ KK) by the QPC 
model to obtain = 0.8 /((tQ^ and rf = 9.6/(crQ^. Com¬ 
paring the above ratios with the experimental data = 3.2 
and rf - 30, they found that ctj a; 1/ V3 is a reasonable es¬ 
timate. In our work we adopt (j)^^ - (iiu + dd + ss)l V3 as 
the flavor S 1/(3) singlet wave function. To show the S 1/(3) 
breaking, the creation strength for mm and dd is defined as 7 , 
while for the strength of ss creation we have 7 / V3, which is 
equivalent to the description in Ref. [17]. In Eq. ( 8 ), sym¬ 
bols p 3 and p 4 stand for momenta of quark and antiquark, 
respectively. J/r,„(p) = lpl^Tf»!(p) is the solid harmonic poly¬ 
nomial is the spin triplet state. The quantum num¬ 

ber of a quark-antiquark pair is = 0^^ determined by 
coupling the orbital angular momentum with the spin angu¬ 
lar momenta, which indicates the conservation of angular mo¬ 
mentum J, P parity and C parity in the course of strong in¬ 
teraction. 0^“^ = (mm + dd + ss)l V3 and , = dij/ V3 are 
the S 1/(3) flavor and color functions, respectively, with i and 
j being the color indices. 

The transition matrix of a process A —> BC can be ex¬ 
pressed as 

{BC\T\A) = 6\pb + , (9) 

where and pc are the momenta of mesons B and C, respec¬ 
tively. |A), \B} and |C) denote mock states [83]. The mock 
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TABLE III: Comparison of experimental data and theoretical results. Here, we also list the )c values for different models. The notation is 
introduced to express mixing states of 'L/, and Here, the results listed in the last column are calculated by the modified GI model with 
yi = 0.03 GeV which gives the least value among several /r’s. 



n “^‘£j 

Experimental values [1] GI model [13] Modified GI model 

D 

1 ‘5o 

1864.84 + 0.07 

1874 

1861 

D‘ 

1^5i 

2010.26 ± 0.07 

2038 

2020 

Z);(2400) 

l^Po 

2318 + 29 

2398 

2365 

£)i(2420) 

IPi 

2421.4 + 0.6 

2467 

2434 

£>,(2430) 

IPi 

2427 ±26 + 25 

2455 

2424 

£>*(2460) 

1 V 2 

2464.3 ± 1.6 

2501 

2468 

0(2550) 

2‘5o 

2539.4 + 4.5 + 6.8 [2] 

2583 

2534 

0*(2600) 

2^51 

2608.7 ± 2.4 ± 2.5 [2] 

2645 

2593 

0(2750) 

1 O 2 

2752.4+ 1.7 + 2.7 [2] 

2845 

2789 

O*(2760) 

1 ^Oj 

2763.3 ± 2.3 ± 2.3 [2] 

2833 

2779 


- 

- 

45677 

5748 


state of a meson A can be defined as 

Ms,M l 

J dpidp26^(pA - Pi - P2) 

X'i'^LM,(Pl’P2)l^l(Pl)HP2)}, ( 10 ) 

where Xs Ms’ spin, flavor and color wave func¬ 
tions of a meson A, respectively. P 2 ) is a spatial 

wave function of meson A, which can be obtained in the mod¬ 
ified GI model. Furthermore, the amplitude can 

be related to the partial wave amplitude via the Jacob- 
Wick formula [84], i.e., 

M^^iA^BC) - V {LQ-JMja\JaMja) 

x{JbMjb; 


meson. For the charmed mesons which are still missing, we 
use the theoretical predictions calculated in the modified GI 
model, the results listed in Table. 11 and/or Fig. 1, as an input. 


A. Instates 

As the I^Pq state, Dg(2400) has been observed by three 
different experiments. However, the experimental masses are 
quite different from each other as shown in Table I. Therefore, 
in this work we take the mass range (2.29 ~ 2.35 GeV) of 
Dq( 2400) to discuss mass dependence of the calculated decay 
width. Here, Dj](2400) only decays into Dtt. In Fig. 2, we 
present the mass dependence of the decay width of Dq( 2400). 
We And that our results are consistent with experimental data 
F = 276 + 21 + 63 [32], F = 240 + 55 + 59 [39], and F = 
273 + 12 + 48 [40]. 

In the following, we study Di(2420) and Di(2430), which 
are mixtures of the I'Pi and PPi states. Di(2420) and 
Di(2430) satisfy the relation 


Therefore, the total decay width can be expressed as 




'4 JX 


\M 


JL,2 


( 11 ) 


After this brief introduction of the QPC model, we perform 
a phenomenological analysis of charmed mesons in the fol¬ 
lowing. When calculating a decay width, we adopt the numer¬ 
ical wave function for a charmed meson obtained in this work 
and the one for charmed-strange meson from Ref. [12]. Addi¬ 
tionally, we still employ the simple harmonic oscillator wave 
function for light mesons such as n and K, where the corre¬ 
sponding p values are taken from Ref. [41]. We need to em¬ 
phasize that the mass is taken from PDG [1] for the observed 


|£>i(2430)> 

|Di(2420)> 


cos 0\p sin 6ip 
-sin 01 /) COS 01 P 


|l')Pi) 

ypi) 


( 12 ) 


where the mixing angle is Oip - -54.7° which is determined 
by the heavy quark limit [41, 85, 86 ]. 

The width of Di(2420) (F = 70 + 21) was Arst measured 
by the ARGUS Collaboration [23]. However, different exper¬ 
iments provided different widths as F = 13 + 6 / 5 *^ MeV [25], 
F = 23/^/^° MeV [27], and F = 21 + 5 + 8 MeV [33]. Here, 
we take the median as the width, i.e., the measurement F = 
21 +5 + 8 MeV from the Belle Collaboration [33]. Di(2430) is 
a broad state, which has the width F = 384/'°^+75 MeV given 
by Belle [32] and F = 266 + 96 MeV by BaBar [38]. Both 
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FIG. 2: Mass dependence of the width of DJ(2400) and comparison 
with the Belle data [32]. The mass range of Z)q( 2400) is 2920 ~ 2350 
MeV. 

Di(2430) and Di(2420) only decay into Dn. In Fig. 3, we 
show the decay widths of Di(2420) and Di(2430) depending 
on the mixing angle O^p, where our results are consistent with 
the experimental data when taking -45.6° < 9\p < -37.2°, 
which is close to 6\p - -54.7° = - arcsin( V2/3) in the heavy 
quark limit [41, 85, 86]. 



FIG. 3: The 9ip dependence of the decay widths of Di(2420) and 
£>,(2430) . 

The ^2(2460) is considered as the 1^P2 state. Its de¬ 
cay channels are Dn, D*n, Dij and DsK, whose theoreti¬ 
cal values of decay widths are shown in Table IV. The total 
width is 14.86 MeV and the branching ratio 21(^2(2460)° —» 
D+7r-)/S(D*(2460)° ^ is 1.96. There are several dif¬ 

ferent experimental widths for ^2(2460), i.e., T = 20+10 + 5 
MeV from the TPS Collaboration[24], T = 20!j°t®2 '^eV 
from the CLEO Collaboration [27] and 45.6 + 4.4 + 6.7 MeV 
from the Belle Collaboration [32]. Our theoretical width 
is consistent with the TPS data [24] and the CLEO data 


TABLE IV: The partial and total widths of charmed mesons 
£>•(2460), £)(2550)/£)2(2580), and £)(2600)/£)}(2650). We mark the 
OZI-forbidden channels by If the channels depend on other pa¬ 
rameters, we label them with and show the widths in Fig. 4. All 
values are in units of MeV. 


Channels 

£>•(2460) 

£>(2550)/£)y(2580) £)(2600)/£)}(2650) 

Dn 

9.17 

- 

□ 

Dj] 

0.02 

- 

□ 

D,,K 

1.0 X lO-'* 

- 

□ 

D*n 

4.67 

67.56 

□ 

D*r] 

- 

- 

□ 

d:k 

- 

- 

□ 

D*(2400)7r 

- 

4.09 

- 

£>;(2460)7r 

- 

- 

□ 

£>i(2420)7r 

- 

- 

□ 

£>i(2430);r 

- 

- 

□ 

Total width 

14.86 

71.65 

- 


[27]. Furthermore, our ratio is in good agreement with ex¬ 
perimental values !B{D^n^)l!B{D*^n^) - 2.3 + 0.8 [27] and 
!B{D^n-)l!B(D*^n-) = 1.9 + 0.5 [32] measured by CLEO and 
Belle, respectively. 

In this subsection, Dq( 2400), Di(2420), Di(2430) and 
^2(2460) which are well established as IP charmed mesons 
have been analyzed, which makes us safely adopt the QPC 
model and the related model parameters to further study the 
following charmed mesons under discussion. 


B. 25 and 1£) states 

The D{2550)/D y(2580) [2, 3] is usually considered as a can¬ 
didate of the 2'5o state. The main decay channels of D(2550) 
are D*n and Dg(2400);7r as shown in Table IV, which can ex¬ 
plain why BaBar and LHCb first observed D(2550)/Dy(2580) 
in the D*n channel. The total width is obtained as 71.65 
MeV which is comparable with the lower bound of the BaBar 
data [2] and is smaller than the LHCb value [3]. Consider¬ 
ing this situation, we also suggest to do more precise mea¬ 
surement of the resonance parameters of D(2550)/Z)y(2580), 
which will be helpful for further testing the 2'5o assignment 
to D(2550)/Dy(2580). 

In the following, we study D*{2600) jD*j{2650) [2, 3] with 
— 1A which is a mixture of the 2^S \ and l^Di states. Here, 
D*(2600) and its orthogonal partner satisfy 


' |D*(2600))' 


cos 6s D sin 6s D 

' 12^51>' 



^ - sin 6sd cos 6sd , 



where the mixing angle 6s d is introduced to describe mixing 
between D{2^S i) and D{\^Di). 
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FIG. 4: The 9sd dependence of the total and partial decay widths and 
the ratio r(£>7r)/r(£>*;r) of £>*(2600). 


FIG. 5: The 0|o dependence of the total and partial decay widths of 
£)(2750)/£)y(2740). Here, the vertical dash line corresponds to the 
mixing angle 9w = -50.8°. 


The 6s D dependence of the total width, partial decay 
widths, and ratio S(D*(2600) —> n~)I!B{D*(2600) —» 

D*^n^) of D*(2600) is shown in Fig. 4, where two tiny partial 
widths r(D*(2600) ^ D*K) and r(D*(2600) ^ D;(2460)7r) 
are not listed. When taking the range -3.6° < Osd < 1-8°, the 
theoretical ratio is consistent with the BaBar measurement of 
Eq. (2). The obtained total width is about 60 MeV which is 
comparable to the experimental data T = 93+6+13 MeV [2] . 
We also find that the main decay modes of D*(2600) is Dn 
(9 ~ 15 MeV) and D*n (32 ~ 38 MeV), which also explains 
why D*(2600) was first reported in these two decay channels 
[2]. We need to stress that the small mixing angle Osd is ex¬ 
pected because of the large mass difference between D*(2600) 
and its partner D*(2760), which is consistent with the sugges¬ 
tion in Ref. [13]. 

As the D-wave charmed meson with = 2 , 

D(2750)/Z)7(2740) is probably either 1D(2") or 1D'(2“) state, 
which satisfies the following relation 


' |1£>(2-)) ' 


J1D'(2-))^ 

v 


cos 6io 
- sin 017 ) 


sin 017 ) 

cos 017) 


| 1 '£» 2 ) 


(14) 


where dio is the mixing angle and can be fixed as dio — 
-50.8° = - arcsin( V3/5) in the heavy quark limit [41, 85, 
87]. 

The decay modes of 1D(2 )/lD'(2“) are shown in Table V. 
If D(2750)/D7(2740) is 1D'(2“), the mixing angle dependence 
of the corresponding partial and total decay widths is given in 
Fig. 5. The range of a mixing angle is obtained as -73.8° < 
6id < -35.7° in Fig. 5 so that the calculated total width is 
consistent with a central value of the experimental data, which 
includes the above heavy quark limit, Oid - -50.8°. Thus, a 
1D'(2“) state is suitable for D(2750)IDj(2740). In addition, 
the main decay modes of D(1D'{2~)) are predicted to be D*n 
(10-25 MeV), Dp (37 ~ 55 MeV), Dw (12 - 17 MeV) and 
D*(2460);r (0 ~ 25 MeV). 


TABLE V: The calculated partial and total decay widths of ID states 
with several possible assignments. If a decay channel is forbidden, 
we mark it by As for the channels which depend on other pa¬ 
rameters, we label them with and show the widths in Figs. 4 and 
5. All values are in units of MeV. 


Channels 

l^Hi 

lD(2-)/lD'(2-) 

I'Ds 

Dn 

76.13 

- 

8.47 

Dri 

9.01 

- 

0.31 

DsK 

11.66 

- 

0.17 

D*n 

35.16 

□ 

7.05 

D*77 

2.68 

□ 

0.11 

d:k 

2.92 

□ 

0.04 

Dp 

26.34 

□ 

0.61 

Dw 

8.87 

□ 

0.21 

D*(2400):/r 

- 

□ 

- 

D*(2460);r 

0.56 

□ 

0.63 

Di(2420)n 

211.72 

□ 

0.21 

Di(2430)n 

0.007 

□ 

0.26 

Total width 

385.06 

- 

18.07 


We need to mention that the widths of 1D'(2“) can be easily 
transformed into those of 1D(2 ), since the width expression 
for 1D'(2 ) with mixing angle dio is equal to that of 1D(2 ) 
with the mixing angle Oid + 90°. Here, we give more pre¬ 
dictions for the missing D{1D(2~)). Its total width can reach 
265 ~ 290 MeV and its main decay modes are D*n (96 -110 
MeV) and ^2(2460)71 (135 - 160 MeV) when taking the 
range -73.8° < 6\d < -35.7°. Here, we take the mass 2737 
MeV for D(1D(2~)) as an input, which is taken from the ex¬ 
perimental data of D 7(2740) [3]. 
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There are two possible assignments of D*{2160)ID*j(2160), 
i.e. either an orthogonal partner of D*(2600)) as given by 
Eq. (13) or a state. Since the mixing angle O^d defined 
in Eq. (13) is quite small, -3.6° < Osd < 1-8°, as shown in 
Eig. 4, so it is legitimate to consider that D*(2760)/D}(2760) 
is dominated by a pure l^Di state. 

The decay modes of D{\^D\) and D{\^D^) states are 
listed in Table V. The theoretical widths of the 
and D{\^D^) states are 385.06 MeV and 18.07 MeV, respec¬ 
tively, both of which deviate from the experimental data, 
E = 60.9 + 5.1 + 3.6 [2] and E = 74.4 + 3.4 + 37.0 MeV 
[3]. This is the obstacle to assign the D{\^Di) and D{\^Di) to 
D*(2760)/D*(2760). 

Apart from the problem of their total decay width, 
in this work we compare the ratio S(D*(2760)° —> 
D+;r-)/S(D(2750)° ^ D*^n-) with the BaBar data [2]. Here, 
Z)(2750) is considered as the 1D'(2“) state. If we assign 
Z)*(2760)/Z)y(2760) to the 1 ^D 3 state, the theoretical ratio is 
0.34 ~ 0.86 with the -73.8° < 6\o < -35.7° range, which is 
in good agreement with experimental measurement [2]. As a 
consequence, D*{2760)ID* {2760) can be tentatively identified 
as the state. However, the future experimental measure¬ 
ment of the resonance parameters for D*(2760)/£)*(2760) is 
an important topic, which will be helpful for giving a final 
and definite answer to the assignment of D*{2760)ID*j{2760). 

Our calculation also shows that the orthogonal partner of 
D*(2600) should be a broad state with the total width 385 
MeV and its main decay channels are Dn, D*n, Dp and 
Di{2420)n. Here, the ratio 


!B{D*'{1-) D*7i) 

S{D*'{\-) -> Dn) 


is also predicted^. 


(15) 


C. 35 states 


We further discuss D}(3000) as D(3^5i). In Table VI, the 
total width of D(3'5 1 ) is 80.36 MeV which is pretty close to 
the lower limit of experimental measurement for D*(3000), 
E = 110.5 + 11.5 MeV [3]. The main decay modes of 
D(3*5'o) are Djt, D*7i and Dp, where the partial decay width 
of D{3^S i) —> Dn contribute about 17% to the total decay 
width which naturally explains why the observed decay mode 
of Dj{3000) is Dn. By the above study, we conclude that 
D*{3000) as D{3^Si) is suitalbe. In addition, we also predict 
the ratio 


!B{D{3^Si) Dn) 
S{D{3^Si) D*n) 


(16) 


which can be tested in future experiments. 

In this work we do not consider the mixing effect between 
2^5 1 and 3^5 1 since the mass gap between these is as large as 
400 MeV as shown in Table 11, which is a different situation 
from the case of mixing between 2^5 1 and l^Di discussed in 
Sect. IV B. 


D. IT states 

There exist four IT states in the charmed meson family. We 
list their decay behaviors in Table VI. 

As for the D{1^F2) state, the mass by the modified GI 
model is 3053 MeV predicted in Table II. The D*j{3000) re¬ 
cently observed can be a possible candidate of D{ 1 ^F 2 ) ac¬ 
cording to the mass spectrum analysis. The main decay chan¬ 
nels of D{1^F2) are Dn, D*n, Dp and Di(2420)7: and the total 
width is 222.02 MeV, which is about two times larger than 
the experimental data, T = 110.5 +11.5 MeV [3]. Therefore, 
D*j{3000) is not a suitable candidate for D{\^F 2 ). Here, the 
ratio 


g(Z)(l^T 2 ) ^ Dn) 
S(D(PT 2 ) ^ D*n) 


(17) 


The mass spectrum analysis suggests the D{3^Sq) assign¬ 
ment for Dy(3000) (see the discussion in Sect. III). Under 
this assignment, we present the decay behaviors of Dj{3000) 
in Table VI. Here, its total width is 90.28 MeV which is 
comparable to the lower limit of the experimental data E = 
188.1 ± 44.8 [3]. The partial decay width of D(3*5o) ^ D*n 
is 43.17 MeV which contributes almost 50% to the total de¬ 
cay width. This is consistent with the fact that Dy(3000) is 
observed in the channel D*n. As a consequence, Dy(3000) 
can be a good candidate of D(3'5'o). Besides D*n, Dp and 
DQ{2400)n are other two important decay modes. 


^ We also notice recent results of LHCb, where a new resonance Z)*(2760) 
[88] with spin-1 was observed. Its resonance parameters are 

M = 2781 ± 18± 11 ±6MeV, T = 177 ± 32 ± 20 ± 7MeV, 

which are comparable with the theoretical mass (2762 MeV) and width 
(385 MeV) of D*'(T). This means that Z)|(2760) can be the orthogonal 
partner of Z)*(2600). 


is also obtained, which provides a crucial information to fur¬ 
ther test whether the D{1^F2) assignment to D*j{3000) is rea¬ 
sonable. 

D(1T(3^)) and D(1T'(3^)) are mixtures of the 1*^3 and 
1 ^T 3 states, i.e.. 


' |1T(3+)) ' 


cos 01 f sin 01 f 

'|1'T3>' 

JIW)), 


-singly COS 0 iy ^ 



where the mixing angle 9if can be fixed as Oif - -49.1° = 
- arcsin(2/ V7) in the heavy quark limit [41, 85] when further 
discussing their decay properties. 

The predicted mass of D(1T(3^)) is 3032 MeV. If we take 
the assignment D(1T(3^)) to Dy(3000), the obtained total de¬ 
cay width is 187.20 MeV, which is in good agreement with the 
experimental measurement E = 188.1 + 44.8 MeV [3]. The 
main decay modes are D*n and ^2(2460)71, which can explain 
why Dy(3000) was first observed by experiment in the D*n 
channel [3]. As a consequence, Dj{3000) can be reasonably 
regarded as the D{1F{3^)) charmed meson. 
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TABLE VI: The calculated partial and total decay widths of Dj(3000) and Dj(3000) with several possible assignments. If a decay channel is 
forbidden, we mark it by All mixing angles are given in the heavy quark limit and all values are in units of MeV. 


Channels 

3'So 

3^51 

Ha 

1F(3+) 

IF'(3+) 

HF4 

2^Po 

2F(1+) 

2F'(1+) 

2^P2 

Dn 

- 

13.53 

26.09 

- 

- 

4.97 

72.51 

- 

- 

1.46 

Dri 

- 

1.37 

2.76 

- 

- 

0.29 

6.96 

- 

- 

0.003 

Di 

- 

0.54 

1.71 

- 

- 

0.02 

1.11 

- 

- 

0.17 

D,K 

- 

2.01 

2.78 

- 

- 

0.15 

7.46 

- 

- 

1.2 X 10-* 

D^n 

43.17 

25.68 

18.83 

42.23 

8.39 

5.31 

- 

61.73 

6.46 

0.12 

D-n 

2.37 

1.78 

1.76 

3.75 

0.35 

0.21 

- 

4.30 

1.26 

0.23 

DW 

0.004 

1.1 X 10-^ 

0.06 

4.9 X 10-* 

8.8 X 10-® 

7.4 X 10-5 

- 

0.90 

0.001 

0.04 

d:k 

3.10 

2.24 

1.53 

3.06 

0.16 

0.09 

- 

3.77 

1.33 

0.26 

Dp 

10.22 

12.16 

19.71 

4.10 

42.86 

2.95 

- 

36.44 

8.80 

13.57 

Dw 

3.26 

4.31 

6.80 

1.28 

12.99 

0.94 

- 

12.56 

2.75 

4.34 

D,K- 

0.004 

0.71 

0.47 

0.01 

0.64 

0.007 

- 

1.94 

6.75 

0.72 

D-p 

2.54 

0.02 

6.41 

16.62 

16.41 

56.94 

138.25 

41.94 

23.78 

5.60 

D'cj 

0.83 

0.006 

2.29 

5.08 

5.45 

18.37 

43.89 

13.95 

7.78 

1.81 

d:k- 

- 

1.1 X 10^^ 

6.3 X 10-^ 

- 

- 

5.8 X 10-* 

3.59 

- 

- 

4.29 

£)*(2400);r 

11.27 

- 

- 

0.23 

0.98 

- 

- 

7.08 

14.01 

- 

D;i2400)ri 

2.79 

- 

- 

0.002 

0.005 

- 

- 

0.39 

0.90 

- 

A, (2317) a: 

5.05 

- 

- 

0.004 

0.02 

- 

- 

0.84 

2.09 

- 

A*(2460);r 

5.67 

6.93 

11.97 

109.45 

3.27 

2.37 

- 

81.44 

9.46 

5.51 

L»i(2420);r 

- 

4.04 

116.51 

1.26 

1.82 

0.33 

11.78 

13.03 

6.85 

7.08 

Ai(2420);; 

- 

0.03 

2.21 

0 

0 

4.0 X 10-* 

0.32 

0.01 

0.008 

0.02 

Di(2430)n 

- 

1.32 

0.31 

0.12 

0.41 

1.55 

11.64 

9.00 

5.03 

21.36 

Dii2430)ri 

- 

1.40 

6.8 X 10-“* 

0 

0 

5.1 X 10* 

0.27 

0.006 

0.003 

0.43 

D,(2460)K 

- 

2.28 

0.002 

3.3 X 10-« 

1.2 X 10“^ 

0.002 

0.57 

0.08 

0.05 

0.88 

Total width 

90.28 

80.36 

222.02 

187.20 

93.76 

94.50 

298.35 

289.41 

97.31 

68.89 


As an orthogonal partner of D(1F{3^)), the theoretical 
mass of D(1F'(3^)) by the modified GI model is 3063 MeV 
which is about 100 MeV above the experimental data of 
Dy(3000). Thus, we also discuss the assignment 

of Dy(3000). The results shown in Table VI indicate that D*n, 
Dp, Du) and D*p are the main decay channels and the total 
width is 93.76 MeV comparable with the lower limit of the 
experimental data [3]. Therefore, we can not fully exclude 
the possibility of that Dj(3000) is a candidate of D(TT'(3^)). 
In order to distinguish possible assignments D(1F(3^)) and 
D(1F'(3^)) of D/(3000), a precise measurement of the total 
and partial widths of Dy(3000) will be a main task in future 
experiments. 

D(PF 4 ) is a possible assignment to D}(3000) according to 
only the mass spectrum analysis since the mass of D( 1 ^F 4 ) 
is calculated as 3037 MeV close to the experimental data of 
Dy(3000). In Table VI, we list the decay channels of D(1^F4). 
Here, D*p and D*(jj are the main decay channels, and the total 
width is 94.50 MeV which is consistent with the experimental 
data r = 110.5 + 11.5 MeV [3]. However, the decay width of 


D{\^F 4 ) —> D*7i is only 4.97 MeV which contributes 5% to 
the total decay width. Thus, it is difficult to explain why the 
observed channel of 77^(3000) is D*ji. By this analysis, we can 
conclude that D*(3000) is not a good candidate of D{l^F 4 ). 
We also give extra information of the typical ratio, i.e., 


S{D{I^F 4 ) ^ Dn) 

\ -- = 0.94. 

S(D(PF4) ^ D*ji) 


(19) 


although these are not main decay modes. 


E. IP states 

In the following, we discuss whether Dy(3000) and 
Dj{3000) can be categorized into the 2P states. 

In Table VI, we list the decay modes of D{2^Pq), where the 
main decay modes are Dn, D*p and D*u) and the total width 
is 298.35 MeV which is far larger than the experimental data 
r = 110.5 +11.5 MeV [3]. D(2^Po) has mass 2856 MeV 
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predicted in the modified GI model which is about 150 MeV 
lower than the experimental mass of D*(3000). Hence, if we 
consider this mass spectrum analysis together with the present 
calculation of the decay behaviors, the D{ 2 ^Pq) assignment to 
Z)y(3000) cannot be supported by our study. 

D y(3000) as a candidate of D{2P( 1 +)) or D{2P'(1 ^)) is con¬ 
sidered here. As mixed states, D{2P(l^)) and D(2P'{1^)) have 
a relation 


F. 2D states 

The mass of D{2^D\) is 3131 MeV predicted through the 
modified GI model in Table II. We show the decay behaviors 
of D{2^Di) in Table VII, which indicates that D{2^Di) is a 
broad state since the obtained total width is 121.75 MeV. Its 
main decay modes contain Dn, D*n and Di(2420);7r. We also 
predict the branching ratio 


' |2P(U)> ' 


cos 02 P sin 02 P 

'|2'Pi)' 

, |2P'(1+)) ^ 


^ - sin 02 P cos 02 P ^ 



( 20 ) 


S(D(2^Di) ^ D’tt) 
S{D{2^Di) -> Dn) 


( 21 ) 


The charmed mesons, D{2D{2 )) and D{2D'{2 )), satisfy 
the following relation 


with a mixing angle 62 P, where we take 02P = Oip - -54.7° 
[41, 85] to list the numerical results of the decay widths of 
D(2P(1+)) or D(2P'(1+)) in Table VI. 


' |2D(2-)) ' 


, |2D'(2-)) ^ 

v 


cos 02D 
- sin 620 


sin 020 
cos 620 


| 2 '£> 2 ) 

\2?D2) 


( 22 ) 


If Dy(3000) is a D(2f’(I+)) state, D*n. Dp, D*p and 
Z)2(2460);7r are its main decay modes, and its total width can 
reach 289.41 MeV comparable with the experimental width 
r = 188.1 + 44.8 MeV [3]. Although the predicted mass of 
D(2P(1^)) is 2853 MeV, which is about 120 MeV lower than 
experimental value of D y(3000), the above results show that 
there is still a possibility of Dj(3000) as D(2P(1^)). 

If Dj(3000) is D{2P'{ I +)), the obtained total width is 97.31 
MeV which is comparable to the lower limit of the experimen¬ 
tal width of D y(3000). In addition, D*p and DQ(2400)7r are the 
main decay modes. However, D*n is not a main decay channel 
since it only contributes 6.6% to the total decay width, where 
D{2P'(1^)) —» D*7 t occurs via a D-wave in the heavy quark 
limit. By remembering that Dy(3000) was first observed in 
its D*n channel, our results do not favor the D{2P'{1^)) as¬ 
signment of Dy(3000), which is also supported by the mass 
spectrum analysis since the theoretical mass of D(2P'(H)) by 
the modified GI model is 2885 MeV far below the experimen¬ 
tal data. 

D*j(3000) is not a good candidate of the D{2^P2) state, 
which is concluded by the mass spectrum analysis and the 
study of its decay behaviors. In the modified GI model, the 
mass of D( 2 ^P 2 ) is predicted to be 2884 MeV, which is in¬ 
consistent with the mass of £>*(3000). Under the D{ 2 ^P 2 ) 
assignment to Dy(3000), the total and partial decay widths 
are presented in Table VI, where Dp and Di(2430);7r are the 
main decay modes and the total width is 68.89 MeV which is 
comparable to the lower limit of experimental measurement 
r = 110.5 +11.5 MeV [3] released by the LHCb Collabora¬ 
tion. However, the process D{2^P2) —» Dn is a subordinate 
decay channel, which just contributes 2% to the total decay 
width. Accordingly, we conclude that D*(3000) is not a good 
candidate of the D{ 2 ^P 2 ) state. 


where 620 is the mixing angle, which in the heavy quark limit 
we can fix as 620 - -50.8° [41, 85, 87]. 

The D{2D(2~)) has the predicted mass 3115 MeV and has 
the broad total decay width, whose value can reach 111.67 
MeV. Additionally, D*n and D2(2460)7r are the main decay 
channels, and the decay width of D{2D{2~)) D*n can con¬ 
tribute more than 40% to the total decay width. 

With the modified GI model, we get the theoretical mass of 
D{2D'{2^)) to be 3148 MeV. In Table VII, we list the decay 
channels of D(2D'(2~)), the total width is 30.17 MeV, and the 
main decay channels contain Dp and Deo. 

The masses of D{2D(2~)) and D(2D'(2 )) are similar to 
each other and have the same decay modes but with difer- 
ent values. However, we can still distinguish them in exper¬ 
iments in two aspects. Firstly, D(2D(2~)) and D(2D'(2^)) 
are broad and narrow states. Secondly, in the heavy quark 
limit D(2D(2 )) —» D*7r is a purely P-wave decay while 
D{2D'{2^)) —> D*n is a purely P-wave decay [87], which 
is the reason why the total decay widths of D(2D(2~)) and 
D{2D'{2^)) are largely different. 

The mass prediction of D{ 2 ^D 2 ,) is 3129 MeV in the mod¬ 
ified GI model. Since the calculated total width of D{ 2 ^D 2 ,) 
is 29.33 MeV, D{2^D^) is a narrow charmed meson, where 
Dn, ^2(2460)71, Di(2420)7r and Di(2430)7r are the main de¬ 
cay channels. We also obtain the ratio 


!B{D{2^D2,) D*n) 
S{D(2^D3) ^ Dn) 


(23) 


which can be tested in future experiments. 


V. SUMMARY 


In Table VI, we predict abundant information of decay be¬ 
haviors of the 2P states in the charmed meson family, which 
provides valuable hints to search for the missing 2P charmed 
mesons and to test these meson assignments to D*(3000) and 
Dy(3000). 


Looking at the the observed charmed mesons shown in Ta¬ 
ble I, the charmed meson family has become more and more 
abundant, which has stimulated us to perform a more system¬ 
atic phenomenological analysis of higher radial and orbital ex¬ 
citations in the charmed meson family with our great interest. 
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TABLE VII: Decay behaviors of four 2D charmed mesons. Values 
are in units of MeV. 


Channels 

2^Di 

2D(2-) 

2D'(2-) 

2^Di 

Dn 

36.10 

- 

- 

3.09 

Dll 

3.49 

- 

- 

0.11 

Dri' 

2.13 

- 

- 

0.003 

D,K 

3.46 

- 

- 

0.06 

D'tt 

13.26 

47.39 

0.29 

0.89 

D'r] 

0.99 

3.86 

0.02 

0.001 

D'ri' 

0.02 

0.28 

0.13 

0.03 

d;k 

0.74 

3.51 

0.04 

1.3 X 10^^ 

Dp 

2.19 

5.36 

14.75 

1.02 

Dcj 

0.71 

1.73 

4.76 

0.32 

D,K- 

0.15 

0.28 

0.01 

0.11 

D'p 

0.01 

0.29 

1.35 

10.11 

D'cj 

0.004 

0.09 

0.42 

3.35 

d;k* 

0.41 

0.47 

0.50 

0.39 

D*(2400);t 

- 

2.54 

2.89 

- 

D*(2400)77 

- 

0.18 

0.31 

- 

D*^(2317)X 

- 

0.24 

0.36 

- 

D;(2460);r 

7.19 

39.92 

0.60 

2.38 

D*2(2573)X 

0.07 

0.37 

0.02 

0.02 

Di(2420);r 

43.32 

1.41 

0.81 

2.24 

Di (2420)77 

0.63 

0.10 

0.12 

0.03 

D,(2430);r 

6.05 

3.36 

2.28 

4.65 

Di (2430)77 

0.27 

0.12 

0.14 

0.25 

Dji(2460)X 

0.34 

0.14 

0.16 

0.27 

Dji(2536)X 

0.22 

0.03 

0.05 

0.007 

Total width 

121.75 

111.67 

30.17 

29.33 


In the present work we have done two major tasks. Firstly, a 
mass spectrum analysis has been given by adopting the mod¬ 
ified GI model, where the screening effect is taken into ac¬ 
count. Secondly, the OZI-allowed two-body strong decays of 
charmed mesons under discussion have been obtained via the 
QPC model. 

In this work, we have revealed the underlying structures 
of the observed charmed states D(2550), D*(2600), Z)(2750), 
D*(2760), Dy(2740), Z)}(2760), Z)y(3000), and D}(3000). In 
Table VIII we give a summary of possible assignments for 
D*j(3000) and Dy(3000). Comparing the theoretical and ex¬ 
perimental data, we rate the possibility of each assignment 
with stars. Additionally, we have provided more abundant 
properties of these particles including some typical decay ra¬ 
tios and partial decay widths, which are critical to test these 
possible assignments of charmed mesons. 

In the following years, exploration of higher radial and or¬ 
bital excitations in the charmed meson family will be one of 
the main projects in Belle, LHCb, and forthcoming Bellell. 
In this work, we have also predicted some missing charmed 
mesons, where their masses and decay behaviors have been 
provided. This information is helpful for experimental study 
of the missing states in the charmed meson family. We also 
expect more experimental observation of charmed mesons in 
future. 
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